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Our numerical result suggests that the dispersion relation of the plasmino mode has a minimum 
at nonzero momentum even near the critical temperature. The quark self-energy is also analyzed 
by using the analyticy of the inverse propagator, which is found to be consistent with the spectral 
function estimated by the two-pole ansatz. 
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1. Introduction 

At asymptotically high temperature (T), properties of strongly interacting matter described by 
Quantum Chromodynamics (QCD) can be calculated using perturbative techniques. In this limit, 
it is known that the collective excitations of quarks develop a mass gap (thermal mass) and a decay 
rate proportional to gT and g 2 T, respectively, where g denotes the gauge coupling [p. Furthermore, 
the dispersion relation splits into two branches, the normal and plasmino modes. As the decay rate 
parametrically grows faster than the thermal mass as g increases, it is naively expected that the 
quark quasi-particles cease to exist as T is lowered. On the other hand, quark number scaling of 
the elliptic flow observed in RHIC experiments indicates the existence of quasi-particles having a 
quark quantum number To understand properties of the matter near T c , especially the quasi- 
particle nature of elementary excitations, therefore, it is desirable to explore the spectral properties 
of quarks within nonperturbative techniques. 

Recently, the correlation function of quarks at nonzero T has been analyzed on the quenched 
lattice with size up to N%xN T = 64 3 x 16 in Landau gauge [Q]. In these studies, it is found that the 
quark correlation function above T c obtained on the lattice is well reproduced by the two-pole ansatz 
for the spectral function, where the two poles correspond to the normal and plasmino modes. In the 
chiral limit these modes have identical quasi-particle masses, which are identified to be the thermal 
mass, that are approximately proportional to T. These calculations, however, also showed that the 
quasi-particle masses, which perturbatively arise through the resummation of infra-red sensitive 
loops [jlj], are strongly dependent on the physical volume, V, used in the lattice calculations. 

The purpose of the present study is to extend the analysis in Refs. ^ to much larger spatial 
volume, xN r = 128 3 x 16 [Q]. This analysis enables us to investigate the spatial volume de- 
pendence in more detail. The large spatial volume also allows to directly analyze the momentum 
dependence of excitation spectra, i.e. the dispersion relations, more precisely. To obtain better 
understanding on the spectral property of quarks, in addition to the analysis of the spectral func- 
tion, we also try to examine the quark self-energy by numerically taking the inverse of the quark 
propagator. 

2. Quark spectral function and fitting ansatz 

Excitation properties of the quark field are encoded in the quark spectral function p^ v (co,p), 
with /i and v denoting Dirac indices. The Dirac structure of p llv {(0,p) at finite temperature is 
decomposed as 

p flv ((o 1 p) = po(co 1 p){f)^ v -p v ((o,p){p-Y) fiV + p s ((o,p)l flv , (2.1) 

where p = \p\ and p = p/p. In the present study we consider the spectral function above T c for 
two cases; (1) at zero momentum, and (2) in the chiral limit. In these cases, the Dirac structure 
of pn V (G),p) are decomposed by using projection operators With p = 0, p v (co,p) vanishes in 
Eq. (^Tj) and p llv (co,p = 0) is decomposed with the projection operators L± = (1 rfc yj/2 as 

p(co,0) = p¥(co)L + f + p M (co)L_f. (2.2) 
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In the chiral limit and for T > T c , the system possesses the chiral symmetry and p s (co,p) vanishes. 
In this case, p [lv ((0,p) is decomposed with the projection operators P±(p) = (1 ± Y°p-f)/2 as 

p((0,p) = pl((0,p)P + (p)f + p p ((0,p)P-(p)r°. (2.3) 

In order to extract p^ v (oo,p) from lattice QCD simulations we have analyzed the quark corre- 
lation function in Euclidean space 

5 MV (T,p) = iy d\d 3 y^< x ^H^(r,x)Y v (0,y)), (2.4) 

on the lattice in quenched approximation in Landau gauge. Nonperturbatively-improved clover 
fermion is used for the analysis. Equation ( |2~4| ) is related to the spectral function as 

/■oo e (i/2-tr)to/r 
S»v(i,p) = J_jG> ea/2T + e _ a/2T Ptiy(<°iP)- ( 2 - 5 ) 

In order to examine the quark spectral function from lattice correlator, we follow the approach 
taken in [gj], which makes use of a two-pole ansatz for the spectrum, 

pf ' P ( to) = Z x 8 (co - E x ) + Z 2 8 ( co + E 2 ) . (2.6) 

Here, Zi 2, and £12 > 0, are fitting parameters that will be determined from correlated fits to the 
lattice correlator: and represent the residues and positions of poles, respectively. When 
comparing the fit results with spectral functions obtained in perturbative calculations one can iden- 
tify the pole at 00 = E\ to be the normal mode, while the one at 00 = —E 2 corresponds to the plas- 
mino mode []5], §]. To determine the fit parameters with correlated fits, we use lattice data points 
at T m i n < xja < N T — T m ; n with T m ; n = 4. We found that the ansatz Eq. ( ^^ ) gives a reasonable 
chi-square, # 2 /dof ~ 1, over wide parameter ranges see, however, Ref. |4j] which discusses 
problems associated with the analysis of the spectral function with an ansatz. 



3. Thermal mass and dispersion relations 

The quark spectral function is analyzed (1) for p = as a fuctions of bare quark mass, mo, 



and (2) for mq = as a fuctions of p. For each case, we use the projection Eqs. (2.2) and (2.3), 
respectively. 

In the analysis of p± I (ft)) at p = on the largest lattice with No/Nt = 8, we found that the 
mo dependence of fitting parameters qualitatively agrees with previous results obtained on lattices 
with N a /N T = 4 and 3 in Ref. [J3J] . From the mo dependence of p±(oo) one can define the critical 
hopping parameter for the chiral limit and the thermal mass, mj [||, Q]. In Fig. [j], we show the value 
of mj with N a /N t = 8, 4 and 3 for T/T c = 1.5 and 3. To infer the thermal mass in the infinite 
volume limit, we performed an extrapolation of mj to infinite volume with an ansatz m^l/V) ~ 
m7-(0)exp(c/V) using the results with N a /N t = 8 and 4. The result of this extrapolation is shown 
in Fig. fil. Within the extrapolation, the value of mj{\/V — > 0) coincides the one obtained with 
N a /N x = 8 within the statistical error. 

Next, we set mo = and analyze the momentum dependence of the excitation spectra at 
nonzero momentum using the decomposition given in Eq. (|2J|). In Fig. we show the momentum 
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Figure 1: Extrapolation of the thermal mass of the quark mj obtained by the pole ansatz to the infinite 
volume limit for T /T c = 3, 1 .5. 




Figure 2: Dependences of the fitting parameters E\ and Ei and the ratio +Z2) on the lattice momen- 

tum p = (1 1 a) sin(pa) for T/T c = 1.5 and 3. See, the text for the explanation of other lines. 



dependence of E\ and E% normalized by mj and Z2/ (Z\ +Z2) for T/T c = 1 .5 and 3. The horizontal 
axis represents the momentum of free Wilson fermions on the lattice, p = (1 /a) sin pa, normalized 
by mj. The figure shows that E\ > £2 is satisfied in accordance with the relation between the nor- 
mal and plasmino dispersions at asymptotically high temperature. The figure also shows that the 
value of £2 for the lowest non-zero momentum, p m i„ = 2%{N Z /N G )T , is significantly lower than 
mj. Provided that the value of £2 in our two-pole ansatz represents the dispersion relation of the 
plasmino mode, this result serves as direct evidence for the existence of the plasmino minimum in 
the non-perturbative analysis. One, however, has to be careful with this interpretation because the 
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Euclidean correlator is insensitive to the spectral function at low energy, |ffl| < T, and analysis of 
the spectrum in this energy range has large uncertainty [Qj. 

4. Quark self-energy 

In order to investigate dynamical properties of the system in lattice simulations, one has to 
take the analytic continuation from a Euclidean correlation function obtained on the lattice to a 
real-time propagator. This analytic continuation, however, is a famous ill-posed problem, because 
one has to infer the real-time propagator, which is a continuous function, from finite and noisy data 



obtained on the Monte Carlo simulations. In previous sections, we have used an ansatz Eq. ( |2.5| ) 
for the spectral function to avoid this difficulty. Although such an analysis would be convenient to 
understand a qualitative structure of the spectrum, details of the spectrum is not accessible. Even 
if one uses the maximum entropy analysis which infers the spectral function without introducing 
an ansatz, the resulting spectral image has uncertainty in the analyses of lattice correlators with 
typical statistics. To make the analytic continuation more robust, therefore, it is desirable to have a 



different formula which relates real- and imaginary-time functions besides Eq. (|2.5|). 

Here, we propose to exploit the inverse propagator for this purpose. The inverse of the retarded 
quark propagator, S R ((0,p), is written as 

[S R (co, P )}- 1 = [S«((D,p)}- 1 -Z R (co,p), (4.1) 

where S R (co,p) and L R (co,p) denote the retarded free-quark propagator and self-energy, respec- 
tively. Let us first derive a formula like Eq. ( |2.5[ ) relating [^(ft),/?)] -1 to a Euclidean function. For 
this purpose we first remark that [^(ft),/?)] -1 is analytic in the upper-half complex-energy plane, 
C + , as well as S R ((0,p). This statement is easily verified by the fact that L R (co,p) is analytic in 
C + by definition. Using this property of [,S R (ft),p)] _1 and Kramers-Kronig relation, by taking a 



similar procedure to derive Eq. (2.5) one arrives at a formula which connects the inverse retarded 



propagator to a Euclidean function 

,°o (l/2-T7>/r r (l/2-T7>/T 

(4.2) 

where [S] -1 (t,/?) is the inverse of the Matsubara propagator S(t,p). Remark, however, that 
Eq. Q4.2[ ) cannot apply to z = 0. In the last equality of Eq. ( p~2| ) we have used the fact that 
[So((0,p)}~ 1 is real 1 , and hence lm[S((0,p)]~ l = -lmL R (co,p). 

The inverse propagator [5] _1 (t,/?) is calculated by inverting the correlation function S(x,p) 
obtained on the lattice. Since S(x,p) is block-diagonal in frequency space, this inversion is most 
conveniently taken after the Fourier transformation, i.e. 

[S]- 1 (r,p) = r£ [S(^,/>)]-V ! ™S (4.3) 

n 

with Matsubara frequency (On = (2n + \)%T, where [ ]~ ! in the r.h.s. represents an inverse of 4 x 4 
matrix for Dirac indices, while that in the l.h.s. means the inverse of the whole propagator. Since 



On the lattice, however, [Sq(cq,p)] takes the imaginary part as a discretization effect. 
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Eq. (4.2) has the same form as Eq. (2.5), once [5] _1 (t,p) is constructed one can infer the real-time 
self-enegy lmL(co,p) with the same techniques to analyze the spectral function with Eq. (2.5), such 
as the maximum entropy method. 

Remarks in constructing [S]~ l (z,p) are in order here. First, [S]~ l (z,p) is not the thermal 
average of the fermion matrix K = ilp — m. This is because the propagator in the l.h.s. of Eq. (1.1) 
is defined by the thermal average of K 1 . The thermal average thus must be taken for the inverse 
of K. Second, to obtain [S]~ l (z,p) one needs all elements of S(z,p) including the value at z = 0. 
For correlators having a positive mass dimension, S(z,p) is singular at z = in the continuum 
limit, and one must exclude this point from the analysis. The quark correlator, on the other hand, 
has zero mass dimension and takes a finite value even at the origin. No difficulty thus in principle 
arises with the use of this point. The use of S(z,p) near the source, however, is troublesome 
because correlators with small z/a receive strong lattice artifacts due to the overlap of operators on 
the lattice. We will later investigate this effect by directly analysing [S]^ 1 (z,p). 

The analysis of [S]~ l (z,p) and lmL R (co,p) has more advantages. First, in the perturbative 
analysis of quark propagator one must first calculate the self-energy. In this sense, the self-energy 
is more useful and fundamental quantity than the propagator in the analytic studies. If one wants 
to compare the analytic result with the lattice one, the comparison is usually made in terms of the 
propagator. The data of the inverse propagator, [S]^ 1 (z,p), on the other hand, enables to make this 
comparison in terms of the self-energy. Next, microscopic physics behind the spectral properties 
would become more apparent through the analysis of ImZ(ft),/?). Because the value of ImE(ft),p) is 
related to elementary processes via the optical theorem, one can give a direct interpretations to the 
structure of lmL(co,p). Of course, p(co,p) and lmL(co,p) should have one-to-one correspondence 
under a given regularization, and in principle one of these functions must contain all information of 
excitation properties. In lattice simulations, however, one cannot determine the real-time function 
definitely. The analysis of a different function thus can make different physics more apparent. 

Because the same correlator as in Eq. (2^5) is used in the construction of [5] _1 (T,p), one 
may think that Eq. {<^2) does not provide any new information on the real-time function besides 
Eq. (2^). We, however, remark that Eq. (O) encodes the analyticy of the inverse propagator in 
C + , which is not taken into account in Eq. (2^). An appropriate use of Eq. (12) thus should enable 
us to extract additional information on the real-time function. 

In Fig. ||, we show the inverse quark propagator [S]q 1 (z,0) = Tr[[5] _1 (T,0)y°] in imaginary 
time in the chiral limit for T = 3T C , with several values of N a and N r , which are constructed 
from the same quark correlators analyzed in the previous sections [0, 0J|. Errorbars are estimated 
by the jackknife analysis. In the figure, we also show the values of [S]q l (z,0) corresponding 
to the spectral function estimated by the two-pole ansatz Eq. (2J5) by the thin-solid lines: The 
spectrum in the chiral limit, p ( ft) , 0) = Z ( <5 ( ft) — ) + <5 ( ft) + ra^ ) ) , corresponds to Im [5 ( ft) , 0) ] q 1 = 
ImTr[I(ft),0)y°] = m 2 T /(TZ)8((o). One finds from the figure that [S]q \z,0) is consistent with the 
prediction of the two-pole ansatz within statistics near zT = 0.5. On the other hand, [§]$ 1 (z,0) 
shows a strong deviation from the constant near the source. Similar results are obtained for nonzero 
mo and p. 

The fact that the values of [SJq^t,/?) coincide with the ones predicted by the pole ansatz 
around zT = 0.5 would support (1) the validity of the pole ansatz for the quark spectrum, and (2) 
the relevance of the evaluation of ImE(ft), p) with Eq. (4.2). The large deviation from the constant 
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Figure 3: The quark self-energy in Euclidean space as a function of imaginary time z for T/T c = 3. 



near the source would be attributed to the distortion effects; in fact, the deviation is more prominent 
on the course lattice. The errorbars and the distortion effect with the present lattice data are too 
large to constrain the form of the spectral function with Eq. (L2). Much finer lattice and higher 
statistics are needed to proceed the analysis of the self-energy with Eq. (L2) further. 

To summarize, in this report we analyzed the spectral properties of quarks on the quenched 
lattice with size x N t = 128 3 x 16. We found that the spatial volume dependence of the quark 
thermal mass defined by the pole ansatz tends to converge on this volume. The dispersion relations 
of the normal and plasmino modes are investigated within the two-pole ansatz. The result indicates 
the existence of a minimum of the plasmino dispersion at nonzero momentum. We also introduced 
attempts to analyze the quark self-energy. 

This proceedings is based on the collabolation done with O. Kaczmarek, F. Karsch, W. Soeld- 
ner, M. Asakawa, and S. Takotani. Numerical simulations of this study have been performed on 
the BlueGene/L at the New York Center for Computational Sciences (NYCCS). 
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